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Abstract 

We discuss new types of second order phase transitions in holography by 
constructing striped black holes in D = 4 with AdS& asymptotics. In 
the context of AdS/CFT, they provide the gravity duals to field theory 
phases in which translational symmetry is spontaneously broken due to 
the formation of current density waves. These black holes are associated to 
three dimensional CFTs at finite temperature and deformed by a uniform 
chemical potential. We numerically solve a non-linear system of PDEs 
in order to construct the black hole geometries and extract some of their 
thermodynamic properties. 



1 Introduction 



The AdS / CFT correspondence provides a powerful tool for studying strongly coupled 
systems. Its possible application to problems relevant to condensed matter physics 
[1] has attracted a lot of attention over the last few years. One aspect is the phase 
structure of CFTs at finite temperature and deformed by a chemical potential. Ther- 
mal phase transitions of CFTs can be studied in the context of holography via phase 
transitions of black holes with AdS asymptotics. A canonical example for the bulk 
dual of the normal phase is the electrically charged AdS-Reissner-Nordstrom (AdS- 
RN) black hole. At low temperatures, these solutions can develop different types of 
instabilities and the system will be described by a new branch of black hole solutions 
in general. 

A prototype example of such a phase transition is provided by bulk theories which 
contain scalars that are minimally coupled to the U(l) gauge field used to deform the 
boundary theory by the chemical potential [2]. In this case, when the normal phase 
becomes unstable against fluctuations of a charged scalar it leads to low temperature 
black holes where the global U(l) symmetry is broken. The gauge field becomes 
massive and we are naturally lead to a holographic description of a superconduc- 
tor/superfluid 013111]. 

Another class of instabilities leads to spontaneous breaking of the Euclidian group 
of symmetries preserved after deforming the theory by a uniform chemical potential. 
Such phenomena are common in strongly coupled condensed matter systems [5] and 
are worth investigating in the context of holography. The first examples of spatially 
modulated unstable modes were realized in D = 5 bulk dimensions in [51 [7] (see also 
[8] for a probe setup). In these five dimensional examples, the unstable modes break 
translations but they still preserve a Bianchi VIIq subgroup giving the dual field 
theory a helical structure. The slices of constant bulk radius are still homogeneous 
and the problem of finding the backreacted black holes was reduced to solving a 
system of ODEs in (7J |9]. The transition was found to be of second order, in general, 
while the spatial modulation persists all the way down to zero temperature and the 
corresponding symmetry breaking ground states were identified in 0HD]. Similar zero 
temperature solutions have appeared in a classification of homogeneous branes-like 
solutions in [TT| H2]. 

Modulated instabilities of electrically charged black holes have also been found in 
D = 4 [T3][]and more recently in [T5] for any D > 4. Unlike the homogeneous D = 5 

1 See also [T3] for a probe brane setup 
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examples, in this case the backreacted geometry is expected to be inhomogeneous 
and the solution of PDEs seems necessary. A first construction of such backreacted 
geometries appeared recently in [16] for a particular model of the classes studied in 
[13]. The authors employed a numerical method developed in [T7] to integrate the 
PDEs resulting from the equations of motion and presented solutions for temperatures 
below T rs 0.9 T c . It was found that for a range of temperatures down to T f» 0.1 T c 
the new inhomogeneous solutions where not thermodynamically preferred over the 
AdS-RN and the system remained dynamically unstable. 

Here we will consider two different models, including the one investigated in [16] , 
and we will solve the system of PDEs by employing the DeTurck method [T8l [19] in 
order to solve Einstein's equations. This method was recently used in [2111 EI] in order 
to study transport properties of backgrounds with explicitly broken translations. For 
both models we are considering, we find a second order transition with the striped 
phase being thermodynamically favored for all T < T c . We present solutions starting 
from T 0.999 T c but it is possible to construct solutions arbitrarily close to the 
critical temperature. 

The remaining sections of this paper are organized as follows. In section [2] we will 
introduce the D = 4 models in which we will be interested and review the spatially 
modulated instabilities of the AdS-RN black hole. In section [3j after describing 
the method we chose to use, we will present the resulting solutions along with their 
thermodynamic properties. We conclude in section [4] with a summary and discussion. 



2 The Einstein-Maxwell-pseudoscalar model 



We consider the D = 4 bulk theory which is described by a metric g^ u , a gauge field 
and a pseudoscalar field ip 



S = J d'xV^g (0*0 



2 ^**-vm 



^j-lJe^)FAF (2.1) 



with F = dA being the field strength. The action (|2.1|) yields the equations of motion 

= 
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d{r* F + $F) = 
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(2.2) 



2 



As discussed in [T3j, when the functions V, r and 9 appearing in the action (2.1) 
admit the expansions 



V = -6 + -m 2 s v 2 + ... 



T 



1 



n 
12 



2^ 



</? + .... (2.3) 



around = 0, the equations of motion (2.2) admit the purely electric AdS-RN black 
hole solution 

ds\ =\ ( — f(z) dt 2 + + dx\ + dx\ 
V f{z) 



A =ii (I 
1 
2 



), f = 

f( z ) =l(i- z ) (-jj, 2 z 3 + 4z 2 + 4z + 4) (2.4) 
with fi being the chemical potential at which the dual field theory is held. In these 



coordinates, the black hole (2.4) has a regular horizon located at z — 1 while the 



AdS± boundary is at z = 0. A simple surface gravity calculation gives the temperature 
T = (12 — yU 2 ) / (87r) showing that the system reaches its zero temperature limit when 
= a/12. The near horizon limit of the extremal solution is simply AdSz x ^ 2 with 
a non-trivial background gauge field 



dsi 



12 



ds 2 (AdS 2 ) + dx\ + dx 



F =V12Vol (AdS 2 ) 



(2.5) 



The perturbative analysis of [13] showed that when the, wavenumber k dependent, 



matrix 



M 



( k2 7s k \ 
24V3A; 24 + A; 2 -c x k 

\ -c x k k 2 + m 2 / 



m 2 = m 2 + n 



(2.6) 



has eigenvalues that violate the AdS 2 Breitenlohner-Freedman (BF) bound m 2 < —3, 



the near horizon geometry (2.5) has tachyonic instabilities. Interestingly, for large 



enough values of the coupling Ci, it was found that the lightest tachyonic modes ap- 
pear at finite values of k. In such a situation, it is natural to expect that there will be 
a critical temperature, at the onset of the instability, at which a new inhomogeneous 
branch of black hole solutions will appear. When these solutions exist at lower tem- 
perature and are also thermodynamically preferred, the system will undergo a phase 
transition. These solutions will break the translational symmetries of the normal 
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phase (2.4) down to a discrete group corresponding to the periodicity 2Ti/k of the 
unstable mode. The mode involves the current and momentum densities of the dual 
field theory and since it is vectorial, parity will be broken at the same time. 

Envisaging a continuous transition, its critical temperature can be determined by 
finding the temperatures T(k) at which the corresponding static modes will appear. 
The details of these static modes were studied in |13j and the profile of the temper- 
ature as a function of the wavenumber T(k) confirmed that for the cases where the 
lightest tachyons existed at finite values of k, the maximum of T(k) was also at k ^ 0. 
The critical temperature of the system is precisely at this maximum T c . 

In this paper we will consider two model choices, the first one given by 

V{<p) = -6 cosh (J^ V Y r= cQsh ^^ » 9 = 7 tanh (Vfy) (2.7) 
which has m 2 s = —4, n = 36 and c\ = 6^7- For values 7 > 7 C = 0.9988... 



there exists a wavenumber interval for which one of the mass matrix (2.6 ) eigenvalues 
violates the AdS2 BF bound m 2 (k) < —3. The model corresponding to the choice 
7 = 1 is of particular importance since it can be derived as a consistent truncation 
of D = 11 supergravity compactified on SE 7 manifolds [22| [23j |23]. As the critical 
value 7 C is approached, the critical temperature T c gets close to zero and a numerical 
treatment becomes less accurate. We will take a larger value for 7 = 1.2 for which we 
give a plot of the static mode curve T(k) in figure [I] From this, we deduce that the 
critical wavenumber is k c ~ 1.33 /i while the critical temperature is T c « 0.0598 fi. 
As we will explain later, the plot in figure [T] was made for the quantization choice of 
the pseudoscalar corresponding to dimension A = 2. 

In order to compare the results of the numerical technique we will use with the 
findings of |16j . the second model we will consider has 



V^) = -6-2 V 2 , r = l, = ^V- (2-8) 

A similar plot with the one in figure HI yields T c w 0.012 p and k c w 0.53/2 In 
the next section we will focus on these two models and construct the corresponding 
broken phase black holes for both of them. 



2 Our value for k c differs from the one of jTHj by a factor of y/2 due to the different scale of x\ on 
the boundary. Notice also a typo for the mass of the pseudoscalar -0 in vl of 16 . 
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Figure 1: We plot the temperatures T(k) for which we find a static modulated mode as 
a function of the wavenumber k for the model choice (2.7) . We find that k c « 1.33 \i 
0.0598 /i. 



and T c « 



3 Striped black holes 
3.1 The setup 

In order to find the backreacted black hole solutions, we need to make a consistent 
ansatz for the metric and the gauge field. The natural requirement is that it should 



capture both the normal phase solution (2.4) as well as the perturbative static modes 
of 



ds\ = — 

z z 



when expanded around (2.4). A suitable ansatz for this purpose is 



A 



z h 



-f(z)Q tt dt 2 , fi _ 
z) a t dt + a 2 dx 2 



dz 2 + Qu (dx\ + z 2 Q zl dz) + Q 22 {dx 2 + (1 - z) Q t2 dt)' 



(3.1) 



with all the nine functions appearing in the ansatz T = {Qu, Qzz, Qu, Qzi, Qt2, a t, a>2, h} 
depending on both z and x\ while f(z) is identical with the one that appears in the 



AdS-RN black hole solution (2.4). Since we wish to deal with the constraints of Ein- 



stein's equations by employing the DeTurck method, we have only partially fixed our 



coordinates in (3.1 ). On the other hand, the gauge field ansatz in (3.1 ) leads to second 



order equations for at and a 2 with no further constrains. This method (THJ 49J and 



an appropriate, slightly less general ansatz than (3.1), was recently used in [201 EI] 
in order to explicitly break translations by introducing a lattice deformation of the 



AdS-RN black hole (2.4). 



Instead of solving Einstein's equations in (2.2), we modify it by shifting its left 



hand side by — V^^) with £ M = g uX (T^ x (g) — T^ x (g)) and g being a reference 



metric which should have the same asymptotic and near horizon structure with the 
black holes we wish to construct [HI [19]. This modification turns Einstein's equations 



to the so called Einstein - DeTurck equations. For a metric ansatz like our (3.1 ) they 
are elliptic and after imposing appropriate boundary conditions they should lead to 
a countable^] set of solutions. We will choose g to be the AdS-RN black hole metric 



which is recovered after setting Q tt = Q zz = Qa = 1 and Q z \ = Q t \ = in (3.1) 
giving a space-like £. Our strategy, as in [20, 21J, will be to check that the solution 
generated by solving these equations will have £ 2 ~ at machine precision. For the 
solutions that we will present in this work we will have £ 2 < 10~ 14 . 

In order to numerically solve the system of the PDEs plus boundary conditions we 
will discretize the coordinates z and x\ on a N z x N x grid. The aim is to transform 
the system of differential equations to an algebraic one. In order to approximate 
the derivatives of our functions at the grid points we will use the pseudospectral 
collocation method. More specifically we will represent our functions as a finite sum 
of Chebyshev polynomials for the z coordinate dependence and as a Fourier sum 
for the periodic X\ coordinate. The method we chose to solve the resulting non- 
linear system of equations is the Newton-Raphson method. Typically, we will be 



taking N x = N z = 35 grid points for the first model (2.7). For the second model 



(2.8), we will be mostly interested in getting accurate results for the thermodynamic 
quantities close to the transition. In that case we will choose N x = 28 and N z = 70 
while checking that the highest Fourier mode we resolve stays small down to the 
temperatures we will construct solutions for. 

We now turn to the question of appropriate boundary conditions. In order to 
avoid the introduction of further deformations to the boundary theory, apart from 
the chemical potential /x, we have to set the sources of all other operators equal to 
zero. In particular, for the pseudoscalar if and the spatial component of the gauge 
field a,2, the general fall off near the boundary z = reads 

V? = zip(i)(xi) + z 2 ip {2 )(xi) H 

a 2 = n j +j(xi)z + 0(z 2 ) (3.2) 

For the pseudoscalar we may choose between two different quantization conditions 
corresponding to a dimension A = 1 or A = 2 boundary theory operator. We will 
choose to work with the A = 2 case in which the function (p^) corresponds to its vev 
while (fm corresponds to its source and we should set cpm = 0. Imposing the absence 



3 Whcn a phase transition happens there are at least two solutions with the same boundary 
conditions. 



6 



of current sources we also require that Hj = 0. Under these conditions, the falloff for 
the remaining fields is 

Q tt = 1 + quix,) z 3 + 0(z 4 ), Q zz = 0{z 4 ), 
Q ii = l + q ii {x 1 )z* + 0{z*), i = l,2, 
Qt2 = qt2(x 1 )z 3 + 0(z 4 ), Q zl = 0(z 2 ) 
a t = l-a{x l )z + 0{z 2 ), (3.3) 

for which the equations of motion also impose 

qtt{xi) + qn{x\) + q^Xx) = 

d^qnixx) = 0. (3.4) 



The conditions of equation (3.4) reflect the fact the energy momentum tensor T, 



of the dual theory, which can be found in appendix |XJ is traceless and source free 
i.e. (T^fj) = and <9 M (T^J) = 0. We are therefore lead to imposing the Dirichlet 
boundary conditions 

Qtt(0,xi) = Q«t(0,zi) = Qn(0,xi) = Q 2 2(0,x 1 ) = a t (0,xi) = 1, 
Q,i(0,X!) = Qa(0,xa) = a 2 (0,x 1 ) = h{0, Xl ) = (3.5) 

on the AdS^ boundary. 

Regularity of the solution at the horizon, which is at z = 1, requires that all our 
functions take finite values and admit an analytic expansion in (z — 1) 

7 = JYl, xi) + d z T{l, Xl )(z -!) + ■■■ (3.6) 



Plugging the near horizon expansion (3.6) in the equations of motion (2.2) and ex- 
panding around z = 1 we find a condition of constancy of the surface gravity on the 
horizon Q Z2 (l,xi) = Qa(l, x i) along with eight relations M. [j r ,d z J-',d Xl J-',d 2 l J-') = 
for our horizon data. These nine equations will serve as our boundary conditions 
on the horizon. It is easy to check that under these conditions the temperature of 
the black hole will simply be given by T = (12 — fi 2 ) /8tt. 

We also note that both our equations of motion and our boundary conditions are 
invariant under translations in x\, as they should since are are interesting in breaking 
this symmetry spontaneously. As a result, if F(z,Xi) is a solution of our boundary 
value problem so is T{z, x\ + b) for constant b. Since we will be searching for periodic 
configurations in x\ it is natural to expect that within a period, there will be points 
for which the d Xl derivatives of our functions will separately vanish. In order to fix 
this mode we will in addition require that (1,0) = 0. In the following section 
we report on the solutions we find to the posed problem. 
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3.2 The solutions 



In order to study the thermodynamics we choose to work in the grand canonical en- 
semble where the chemical potential fx is held fixed. In practice, we will be generating 



solutions by varying fx in (3.1) and by using the scaling symmetry of the boundary 
theory we will evaluate all the quantities in units of fx. As previously mentioned, we 
will be focusing on a particular periodicity in x\ for our solutions corresponding to 
the wavelength 2ir/k of the critical mode with k = fc c ~ 1-33 fx for e.g. the first model 



(2.7). For temperatures higher than the critical T c j fx ~ 0.0598, the only solution we 



find to our equations is the RN-black hole (2.4). In accordance with our expecta- 



tions, for temperatures T < T c « 0.0598 fx, we find additional solutions for which the 



functions in our ansatz (3.1) have the anticipated periodic profile in x±. 

For two such solution, with T = 0.6 T c and T = 0.17T C , we present the profiles 
of the functions h and a t in figure [2] for a single period in x\. In [13] it was shown, 
using a perturbative series argument, that the electric component of the gauge field 
a t would be modulated with a period equal to ir/k close to the transition. For both 



models (2.7) and (2.8) we chose to study, a similar argument can be made to all 
orders in perturbation theory in agreement with the periodicity of the right plots in 
figure [2] We also observe that translation breaking effects become stronger in the IR 
as we lower the temperature. 

The averaged entropy and charge densities of our system can be found through 

r2n/k 

S = k ^Q 11 (l,x 1 )Q 22 {l,x 1 )dx 1 
Jo 

kfx f 2n ^ k 

Q = fx H / a(xi) dxi (3.7) 

while the averaged free energy density is 

tt = M - fxQ -T S 

M = 2+^-^ / q tt (xi)dx 1 (3.8) 

with M being the averaged mass density of the black hole, as we explain in Appendix 
|A) In figure [3] we give a plot of the averaged entropy density difference between the 
striped black hole and the RN black hole solution as a function of the temperature. 
We do the same for the averaged free energy density difference in figure [3j We see 
that the new phase of black holes leads to lower entropy as well as free energy. More 
specifically, we find that near the critical temperature A^^|j oc — — and 
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Figure 2: On the left we plot the function h parametrizing the pseudo scalar. On the 
right we plot the function a t . The top plots are for T = 0.6 T c , the bottom plots are 
for T = 0.17T C and have N x = 40, N z = 50. 



A 



oc 



1, compatible with a second order phase transition and the first law 



of thermodynamic^] 



5Q 



-S5T-Q5fi 



(3.9) 



for 5[i = 0. 

Another set of quantities we will now consider is the Fourier modes of the pseu- 
doscalar operator vev 



dn) 
\2) 



k 

2^ 



2-w/k 



-inkx 



1 4>(2)(x 1 ) dxi. 



(3.10) 



Close to the critical temperature T c , the mode with n = 1 is expected to have the mean 



is compatible with 
under consideration 



x(l) 

field behavior -^2. ^ (i _ ^-)V 2 an( i the left plot shown in figure 
such a power law behavior. In more detail, for the specific mode 
with even V(ip), r(ip) and odd 0((p), only modes with odd n will be involved in a 
perturbative expansion of the pseudoscalar. With the n-th mode being sourced at 
n-th order in perturbation theory we expect that we should have the dependence 



'We keep the period of the solutions in the x\ direction fixed. 
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Figure 3: On the left we plot we see that the striped black holes have lower free energy 
density below the critical temperature T c . The right plot shows that the broken phase 



black holes lead to lower entropy. Both plots are for the model (2.7) 
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Figure 4: The left plot shows the first non-trivial Fourier mode of the pseudo scalar 
vev as a function of temperature. The right log-log plot shows its first three non- 
trivial Fourier modes as a function of 1 — T/T c along with a linear fit. 



oc (1 - f ) n ' 2 . In figure 



we present a log-log plot of 



dn) 
\2) 



as a function of 1 — |r 



for the first three non-trivial modes corresponding to n = 1, 3, 5. A linear fit 



log 



\2) 



-oc n + 7„ In 1 



T 

Tc 



(3.11) 



based on our near transition data determines the slopes 71 ~ 0.501, 73 ~ 1.508 and 
75 rs 2.505, in good agreement with our previous argument. Similar behavior is found 
for the momentum ^2(^1) and current density waves. 

We will now examine the charge density wave which forms in our striped phase 
and whose density profile can be read off from a{x\) in (3.3). In figure [5] we plot the 
difference of the averaged charge density AQ from the normal phase as a function of 
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Figure 5: On the left we plot the difference of the average charge density between the 
striped and the normal phase black hole as a function of temperature. On the right 
we show the first non-trivial Fourier mode (n = 2) of the charge density wave. 



temperature. In the same figure we also plot its Fourier mode 

nh f 2n / k 

(n) _ I -mkxi 



q^ = ^- / e~ mkxi a(x 1 )dx 1 (3.12) 
2tt Jo 

with n = 2 as a function of temperature. For both cases we observe the expected 
linear dependence q n=0 > 2 oc 1 — ^ close to the transition as the electric component of 

J- c 

the gauge field a t backreacts at second order. 



We now turn our attention to the model (2.8) and repeat the same procedure 
for k c ~ 0.53 \x and T c f» 0.012 fi. We don't find any qualitative difference between 
the two models close to the transition, as we would expect. The transition is again 
second order and in contrast with the results of [TB], we find solutions with T < T c 
which have lower free energy and entropy than the AdS-RN black hole for a given 
temperature. We plot the differences of the free energy and entropy densities between 
the two phases in figure [6j 

In order to calculate the free energy (3.7) we need to extract the function qu{x\) 
from the falloff (3.3) of the function Q u that parametrizes our ansatz (3.1) e.g. by 
taking a third order derivative at z = 0. Even though the pseudospectral technique 
we have chosen makes the calculation of derivatives of functions quite accurate, it 



still is less accurate than the calculation of e.g. the entropy (3.7) which is using 
directly the values of the functions. The fact that our data satisfies the first law of 
thermodynamics at a good accuracy is quite reassuring. In Appendix [B] we outline 
an equivalent way to phrase the boundary value problem in which the calculation 
of the normalizable modes of our functions is more direct. We have also solved our 
PDEs in the way described there and reproduced the plots of figure |6l 
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Figure 6: We find similar behavior for the striped black holes of the model (2.8). The 
striped phase has lower free energy density and entropy below the critical temperature 
T 

->- c- 



4 Discussion 

We numerically constructed four dimensional black hole geometries with inhomoge- 
neous horizons by backreacting on the instabilities of [13]. These geometries are the 
bulk duals of CFT phases held at finite temperature and deformed by a uniform 
chemical potential. The numerical method we chose to follow [TS1 |T5] gives results 
compatible with a second order phase transition. We only reported on the choices 



(2.7) and (2.8) but since the order of the transition is decided at third order in per- 
turbation theory [22] we expect the character of the transition to remain the same 
for a general class of models. In this work we focused on a particular periodicity 
for the modulation of our solutions but it is clear from the five dimensional homo- 
geneous analogs [9l |26] that the configurations of minimum free energy will have a 
temperature dependent periodicity. The order of the transition is unlikely to change 
due to this phenomenon in general. An open question is that of the minimum free 
energy configuration, it is possible that the thermodynamically preferred mode is a 
superposition in A;— space. To answer this question one would have to solve PDEs in 
three dimensions or follow a perturbative approach fI7\ [2"8~] . 

It would also be interesting to explore the phase transition in theories where the 
near horizon limit of the normal phase black holes at zero temperature is hyperscaling 
violating instead of AdS2 xl 2 . It was recently shown, in [21], that the near horizon 
limit of such black holes can have instabilities of the same nature. Based on the 
general arguments of [30] the phase transition in those theories is expected to exhibit 
qualitative differences reflecting the fact that the normal phase black hole branch has 
zero entropy at its zero temperature limit. 
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An important question is that of the zero temperature ground state. As we 
lower the temperature and the non-linearities of the functions V, r and 9 become 
important, we expect a rich class of ground states to emerge in the IR. In the absence 



of modulation and for a purely electric ansatz for the gauge field, the models (2.1) 
reduce to Einstein-Maxwell-Dilaton. In that thorough classification of possible 

zero temperature IR behavior has been carried out in j3TJ E2J S3]- It would certainly 
be interesting to carry out as a similar classification in an inhomogeneous setting. 
Based on the numerics in [IB] , it was argued by the authors that new IR behavior 
emerges also in D = 4. We plan to further explore the low temperature behavior of 
the striped phases, based on the techniques used in this paper, in the future. 
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A Holographic renormalization 

In order to have a well defined variational problem and a finite total action, we need 



to supply the action (2.1) with appropriate boundary counterterms 



S T = S + J gM ( /C_2V ^"^^) (A,1) 

where the last term ensures that our pseudoscalar has dimension A = 2. We also 
have that JC = lim r ^o 9^ V^n^ with n M an outward pointing unit normal vector and 
7 is the determinant of the induced metric on the boundary at r — > 0. Following 



and given our asymptotic expansion (3.3) along with the condition (3.4), we obtain 
the expression 

(Tu) = V2^2+^- -SquixSj , (T XiXi ) = (l + ^ +3q u {x 1 )j 
(Tt, 2 ) = 4|Sa(zi) (A.2) 

for the non-trivial component of the boundary stress energy tensor. 
The metric of the boundary theory is 

^ s 3 = l{3)fMvdx' 1 dx L ' = —2 dt 2 + dx\ + dx\. (A. 3) 
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After defining the timelike unitary vector u = 2 1//2 d t , the mass and momentum 
densities read 

u 2 

m =^m)U»u u (T^> = 2 + ^ - 3g tt (xi) 



3 

Px 2 =u' I (T llX2 ) = -q t2 (xi). (A.4) 



B An equivalent boundary value problem 



We will now take advantage of the asymptotic falloff (3.3) to make a better choice of 



our functions in order to extract the asymptotic data. Instead of using the functions 



we introduced in the ansatz (3.1), we make the change of variables 



Q mm 0, xi) = 1 + z 3 q mm (z, x x ) , m = t, z, 1, 2 
Q zl (z, xi) = z 2 q zl (z, xi) , Q t 2 (z, %i) = z 3 q t2 {z, Xx) 
a t (z, Xi) = 1 — z a (z, Xi) , a 2 (z, x\) — zj (z, x±) 

If = Z 2 0(2) {z, xi) . (B.l) 

The aim is to solve the problem in such a way that the values of the functions we 
are solving for are directly related to the values of the vevs we wish to calculate. 
Our requirement for a regular horizon at z — 1 leads us again to imposing very 
similar boundary conditions with the ones described in section [3] for the functions 
J~ = |<Zmm> Qrxn <jtx 2 , 3i V?(2) j- F° r example we find that we should impose the 
zeroth law of thermodynamics guaranteed by q u (l,^i) = <jrr (l,^i)- 

The difference in the nature of the boundary conditions we need to impose comes 
from the AdS± boundary at z — 0. Plugging the ansatz in the equations of motion 



(2.2) we find the simple first order relations at z = 



d z q a = -0\ 2 ) + a (o + 2) , d z q rr = ^ j 2 + 3 0\ 2) - ^/i 2 a (a + 2) 
d z qn = -0 2 ( 2 ), d z q 22 = -]-j 2 - 2 {2) , d z q zl = 0( 2 ) d Xl 0(2) 



d z q t 2 = Qt2 + (a + 1) , d z d = a, d z j = 0, d z 0( 2) = (B.2) 

which we can use as our boundary conditions. We also find that that we should 
satisfy the conditions 

q tt (0, xx) + q XlXl (0, Xi) + q X2X2 (0, Xi) = 0, d Xl q XlXl (0, x x ) = 0, q zz (0, x\) = 

(B.3) 
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Given the above boundary conditions for the functions J 7 , we can reproduce the plots 
of figure [6] at relatively low grid resolutions e.g. N x = 25 and N z = 45. We also check 
that the resulting solutions satisfy the constraints (B.3). 
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